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Abstract 

The spectrum of the Dirac operator near zero virtuality obtained in lattice gauge 
simulations is known to be universally described by chiral random matrix theory. We 
address the question of the maximum energy for which this universality persists. For 
this purpose, we analyze large ensembles of complete spectra of the Euclidean Dirac 
operator for staggered fermions. We calculate the disconnected scalar susceptibility 
and the microscopic number variance for the chiral symplectic ensemble of random 
matrices and compare the results with lattice Dirac spectra for quenched SU(2). 
The crossover to a non-universal regime is clearly identified and found to scale 
with the square of the linear lattice size and with f%, in agreement with theoretical 
expectations. 



Recently, it has been shown by several authors that chiral random matrix 
theory (chRMT) is able to reproduce quantitatively spectral properties of the 
Dirac operator obtained from QCD lattice data. This statement is valid both 
for fluctuation properties in the bulk of the spectrum and for microscopic 
spectral properties near zero virtuality, see the reviews [1,2] and Refs. [3-6]. 
This result implies that the spectral fluctuation properties of the Dirac op- 
erator are universal, i.e., determined solely by the underlying symmetry of 
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the problem and quite independent of specific aspects of QCD. The success 
of chRMT poses the question: Which QCD energy scale limits this universal 
behavior? In mesoscopic physics, the analogous scale (i.e., the "Thouless en- 
ergy") is given by E c ~ L~ 2 , where L is the length of the sample. Spectral 
fluctuation properties of a mesoscopic probe obey random matrix theory only 
in energy intervals smaller than Ec- 

Two recent publications [7,8] address the existence of such a scale, here de- 
noted by Armt, in QCD. (Earlier qualitative discussions of the transport prop- 
erties of light quarks in the QCD vacuum can be found in Refs. [9,10], and a 
more quantitative approach was taken recently in Ref. [11].) The scale Armt 
is important since on smaller scales, QCD calculations do not contain system- 
specific information. The authors of Ref. [7] used general arguments and simple 
estimates for Armt, while Ref. [8] provides semi-quantitative results for Armt 
based on the instanton liquid model. It is the purpose of this Letter to deduce 
for the first time values for Armt directly from microscopic QCD lattice data 
and to establish the scaling properties of this quantity both with respect to 
lattice size and coupling constant. A recent analysis of spectral data in the 
bulk [12] yields results which are consistent with our findings. 

We recall that chRMT uses a generating functional of the form 

N f 

Z§ = f D[W\ J] det(D + m f )e-^ tTv( - w ' w) (1) 



with 



V 




(2) 



and a potential, v, which determines the distribution of the matrix elements of 
W. The universal spectral fluctuation properties do not depend on the choice 
of v [13] which is taken to be a Gaussian for convenience, 

v(W ] W) = Y?W ] W . (3) 



In Eq. (1), we consider only the sector of topological charge zero because our 
lattice data agree with the chRMT results in this sector [5]. With N the di- 
mension of the matrix in Eq. (2), E is the absolute value of the chiral conden- 
sate, (ipip) (per flavor). Various gauge theories have different symmetries and, 
hence, different values for (3d. (The index D for Dyson serves to distinguish 
the symmetry parameter from the square of the inverse coupling constant de- 
noted by (3.) For SU(iV c ) and N c > 3 one has /3 D = 2 (chiral Gaussian Unitary 
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Ensemble, chGUE); for N c = 2 and staggered fermions (this is our case) one 
has /3d = 4 (chiral Gaussian Symplectic Ensemble, chGSE); and for N c = 2 
and fermions in the fundamental representation one has (3d = 1 (chiral Gaus- 
sian Orthogonal Ensemble, chGOE), see Ref. [14]. When we apply chRMT to 
quenched lattice calculations, the determinant in Eq. (1) is absent. 

Earlier comparisons have shown that all predictions of chRMT such as sum 
rules, microscopic spectral distributions, spectral correlations in the bulk, 
nearest-neighbor spacing distributions, etc. agree very well with lattice data. 
The single parameter of the model, S, can be determined from the lattice 
data [5] via the Banks-Casher relation [9]. Then, the chRMT predictions are 
parameter free. 

In Refs. [7,8], it was argued that Armt can be estimated with the help of the 
Gell-Mann-Oakes-Renner relation, which yields 

fl 

Armt ~ Y^Tp. ' ( 4 ) 



where f n = 93 MeV is the pion decay constant and V = L 4 is the space-time 
volume. On the lattice, V = Na 4 , where N is the number of lattice sites and a 
is the lattice constant which we set to unity unless otherwise indicated. Using 
the mean level spacing at zero, A = 7r/(EV), Eq. (4) can be expressed in 
dimensionless form, 

A R mt/A ~ -flL 2 . (5) 



To determine Armt and to test the expected dependence of Armt on L and S, 
we use the disconnected spectral susceptibility x dlsc and the S 2 (0, S) statistic. 
In order to avoid confusion between the latter quantity and the value, S, of 
the chiral condensate, we will display the arguments (0,S) in S 2 (0, S). We 
denote the limiting scale for chRMT determined from x dlsc by Armt and that 
determined from S 2 (0, S) by Srmt- We shall see that Armt/ A and 5rmt agree 
within the accuracy of our analysis, although the errors associated with Srmt 
are larger than those for Armt/ A. We shall not address the question of how 
these quantities are related to an intrinsically defined energy scale, cf. Ref. [15]. 

The disconnected spectral susceptibility, x dlsc , is defined in terms of the Dirac 
eigenvalues, A^, obtained in lattice simulations by 



(6) 




N \k~hli (^Afc + m)(i\i + m) 




where the average is over independent gauge field configurations and where m 
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denotes the valence quark mass. Note that for SU(2) all eigenvalues are two- 
fold degenerate; for, e.g., SU(3) the sums would run up to 3N. We study x dlsc 
at zero temperature. Lattice QCD studies of the disconnected and connected 
susceptibilities at finite temperature do exist [16]. In this case, chRMT must 
be supplemented by non-random terms which are model-dependent [17]. Nev- 
ertheless, the universality of the random-matrix results is expected to persist 
for energies below Armt if the (model-dependent) temperature-dependence 
of E is taken into account [18]. However, we shall not address the question of 
finite temperature in this work. 

The sums in Eq. (6) can be written as integrals involving the microscopic 
spectral densities of the Dirac operator, i.e., the spectral densities on the scale 
of the mean level spacing near zero. We have 



X 



disc 



= Au 1 



dx 



Pl (x) 



(x 2 + u 2 ) 2 



dx I dy 



(x 2 + u 2 )(y 2 + u 2 ) 



(7) 



where u = miVE and x dlsc has been rescaled by l/(iVE 2 ) so that all quantities 
in Eq. (7) are dimensionless. The function T2(x,y) is the connected part of 
the microscopic spectral two-point function, p 2 (x,y) = pi{x)p 1 {y) — T 2 (x,y). 
Equation (7) is universal in the sense that all reference to the parameter E, 
which depends on the simulation parameter /3 = 4/g 2 , has been eliminated. 
We now make the transition to chRMT by substituting the random-matrix 
results for the microscopic spectral one- and two-point functions appearing 
in Eq. (7). For the quenched chGSE, we have [19] 

i i 

Pl (x) = 2x 2 j dss 2 J dt [J {2stx)J l {2sx) - tJ {2sx)J 1 {2stx)} (8) 



and 



r 2 (x, y) = (2xy) 2 [S{x, y)S(y, x) + I(x, y)D(x, y)\ (9) 



with 



S(x,y)= J dss 2 J dt[J (2stx)J 1 (2sy) - tJ (2sx)J 1 (2sty)] (10) 



1 1 

I(x,y) = J dss J dt [J (2stx)J (2sy) - J (2sx)J (2sty)] (11) 
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1 1 



D(x,y) = J ds s 3 J dtt[J 1 (2stx)J 1 (2sy) - J 1 (2sx)J 1 (2sty)} , (12) 



where J denotes the Bessel function. A tedious calculation leads to the result 



X 



l 

disc 



4u 2 J ds s 2 K {2su) J dtl {2stu){s{l-t 2 ) 



+AK (2u) [I (2su) + tl (2stu)] - 8stI (2stu)K (2su)} 



Au 2 K 2 (2u) 



J ds I (2su) 



(13) 



where / and K are modified Bessel functions. To the best of our knowledge, 
Eq. (13) presents a novel result. The disconnected susceptibility is the result 
of strong cancellations between the two terms in Eq. (7). For this reason, x dlsc 
is particularly sensitive to deviations from chRMT and well-suited for the 
determination of Armt- 

We turn to a comparison of x dlsc a s predicted from Eq. (13) with lattice data. 
As mentioned above, this is the first time that such a comparison has been 
made. In Ref. [8], Osborn and Verbaarschot presented calculations for x dlsc 
from the instanton liquid model, an effective model for QCD. Their results 
show certain features which are difficult to interpret and which may be due 
to finite-size effects as they suggest. 

Figure 1 shows the dependence of x dlsc on the scaled valence quark mass u, 
defined below Eq. (7), for a typical example, L — 10 and (3 = 2.0. Here and 
below, the values for S are taken from Ref. [5]. (Note that the eigenvalues 
in [5] were measured in units of l/(2a).) The results shown in Fig. 1 were 
obtained without spectral unfolding. We have also unfolded the lattice data, 
but the resulting differences in x dlsc are negligible since the sums in Eq. (6) 
are dominated by small eigenvalues for which the spectral density is approx- 
imately constant. Hence, details of the unfolding procedure are irrelevant for 
the present investigation. 

We note that the uncertainties in the Monte-Carlo data are correlated: The 
entire set of dots in Fig. 1 would shift up or down within the range indicated 
by the error bars if, e.g., the lowest eigenvalue were allowed to move within its 
statistical error. Our interest is focussed on the systematic deviations visible 
above u ~ 7. In order to determine these deviations, we show in Fig. 2 the 
ratio 

ratio = ( X f- ce - X ^ T ) I ( X di - T ) . (14) 
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Fig. 1. The scaled disconnected susceptibility plotted versus the scaled valence quark 
mass. The open squares are lattice data; the dots are the chRMT prediction. The 
data consist of 1416 complete spectra on an N = 10 4 lattice with (3 = 2.00 and 
E = 0.1247. 
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Fig. 2. The relative difference, Eq. (14), of the scaled disconnected susceptibilities 
for the lattice simulation (using the data from Fig. 1 with ./V = 10 4 , (3 = 2.00, and 
S = 0.1247) and chRMT. 

Deviations of this ratio from zero determine Armt/A. The errors in Fig. 2 
are jackknife estimates. Two features in the figure are striking, (i) Below the 
lowest eigenvalue of the Monte-Carlo sample, the errors are too small, (ii) For 
very small values of u, one observes a systematic deviation between the lattice 
results and the chRMT prediction. These features are artefacts of limited 
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Fig. 3. The relative difference of the scaled disconnected susceptibilities plotted 
versus u/L 2 for (3 = 2.00 and four different lattice sizes, iV = 4 4 , 6 4 , 8 4 , and 10 4 . 

statistics and have the following cause. The asymptotic chRMT result for 
very small values of u is 



X 



disc 



(2uf 



--(ln W + 7 )-- 



(15) 



where 7 is Euler's constant. The logarithmic term is generated by the small but 
finite eigenvalue density at small u, see Eq. (7). However, in a given Monte- 
Carlo simulation there is always one smallest eigenvalue, A m ; n . For values of u 
smaller than A m j n iV£, the logarithmic contribution can no longer be obtained 
from the lattice data, see Eq. (6). 

Let mrmt = ti"Armt/A be the value of u at which the strong deviation observed 
in Fig. 2 sets in. According to Eq. (5), wrmt should scale with L 2 [7,8]. To 
check this prediction, we have plotted in Fig. 3 the ratios defined in Eq. (14) 
for L = 4, 6, 8, and 10 as a function of u/L 2 . Obviously, all data fall on the 
same curve confirming our expectation. 

In order to compare results for different values of /3, we note that Wrmt is 
dimensionless but should be proportional to L 2 . The latter quantity should 
scale with a 2 , where a depends on (3. Furthermore, in the scaling regime one 
would expect that E scales with a~ 3 . This suggests that mrmt/ (S 2 / 3 L 2 ) should 
be independent of /3 in the scaling regime. Figure 4 demonstrates that this 
expectation is not supported by the data. It is perhaps not too surprising that 
simple scaling does not work, because the dynamics on the lattice changes in 
a highly complicated manner between (3 = 2.0 and (3 = 2.4. The theoretical 
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Fig. 4. The relative difference of the scaled disconnected susceptibilities plotted 
versus u /(E 2 / 3 L 2 ) for the data of Fig. 3 and additional data for (3 = 2.2, E = 0.0556 
on 6 4 and 8 4 lattices and for j3 = 2.4, E = 0.00863 on a 16 4 lattice. 
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Fig. 5. The data of Fig. 4 plotted versus u/(EL 2 ). 
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Fig. 6. Comparison of the number variance, X 2 (0, S), predicted by chRMT with the 
results for the simulations used in Fig. 3. 

expectation of Eq. (5) is that wrmt should scale with f 2 L 2 [7,8]. A careful 
check of this expectation would require the determination of f n for the lattice 
sizes and (3 values we have used. We have not done this. Instead, we make use 
of the observation [20] that f 2 (in lattice units) scales approximately like S 
for the range of (3 considered here. The plots in Fig. 5 showing the results for 
different (3 versus m/(SL 2 ) support this view. 

Billoire et al. [20] suggest that f 2 = S/3.4 in lattice units. If one interprets 
Fig. 5 as indicating that mrmt/(SL 2 ) is roughly 0.5, this implies (taking into 
account a factor of 1/2 from our normalization of the eigenvalues) that 

A R mt/A « 0.3 f 2 L 2 . (16) 




This result in quenched SU(2) is in agreement with the order of magnitude 
estimate Armt/A ~ f 2 L 2 /n from Refs. [7,8]. 

We now turn to the number variance which is defined as S 2 (0, S) = ((N(I) — 
(N(I))) 2 ) [6]. Here, / is the interval / = [0, S], N(I) is the number of eigen- 
values in J, and the angular brackets denote the ensemble average. In contrast 
to x dlsc , unfolding is important for the T* 2 (0,S) statistic since it leads to a 
significant extension of the length of the interval / for which S 2 (0, S) can be 
determined. We unfolded the spectrum by fitting the unfolding function to 
the average of the spectrum over all configurations. Figure 6 shows that the 
critical value, Srmt, for which deviations from chRMT are observed increases 
with L. In Fig. 7 we see that Srmt decreases with increasing (3 as expected. 
Detailed analysis of all available data sets shows that 
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Fig. 7. Same as Fig. 6 but keeping L fixed and varying f3. 



Srmt ~ (0.3 - 0.7) SL 2 



(17) 



is consistent with the data. Hence, Srmt and mrmt are perfectly consistent. 



The use of the E 2 (0, S) statistic may be conceptually more appealing because 
the analogous quantity in mesoscopic systems is directly related to the Thou- 
less energy. However, our analysis shows that the susceptibility appears to be 
better suited for a quantitative determination of the cross-over point from 
universal to non-universal behavior. 



In conclusion, we have provided the first direct determination of the scale, 
Armt, which limits the validity of random matrix descriptions of lattice QCD. 
This quantity has the correct L 2 -scaling. Moreover, Armt seems to scale rough- 
ly with / 2 as expected on the basis of the Gell-Mann-Oakes-Renner relation. 

It would be very interesting to perform a detailed analysis of the lattice data 
in the diffusive regime, i.e., above Armt, to check the predictions of Ref. [7] 
for this regime and to investigate possible differences between the numerical 
results of Ref. [8] for the instanton liquid model and the lattice data. Such an 
analysis will be the subject of future work. 
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